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6=SEP=19 MTHRTL.SRCJUVXGSINCO.MAR; 1 a1) | 

TITLE UVXSGSINCOS ; Floating Point Sine, Cosine and Sincos | 
; UNCTIONS 

IDENT /2=004/ : File: MTHGSINCOS.MAR EDIT: JCW2004 | 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


1S SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
LY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
Siier OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
HE 
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ES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
TAAHIStERAEO NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
ooRPOhat ite NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONS . 


IBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH 


L 
S NOT SUPPLIED BY DIGITAL. 


Se Se Be Be Oe Oe Oe He Oe Oe Oe Oe Oe Oe He Se Oe Oe Se ee 
ae eeeeeeeeeeee ee eee ene 


® 
* 
® 
® 
& 
© 
® 
® 
® 
* 
® 
® 
® 
® 
® 
* 
® 
® 
® 
* 
® 


LARARABABAASARALEAAAAAAARE SELES EE ESSERE RRA E SELES RRR ERS R RRR R RRR R RRR REESE SD 


FACILITY: MATH LIBRARY 
ABSTRACT: 


MTHSGSIN and MTHS$GCOS are functions which return the pod i, point 
sine or cosine value of their single precision floating point argu- 
ment (radians). The call is standard call-by-reference. 
MTHSGSIN_R7 and MTH$GCOS_R7 are special routines which are the same 
as MTHSGSIN and MTH$GCOS” except a faster non-standard JSB call is 
used with the argument in RO and no registers are saved. 


MTH$GSINCOS is a routine which returns the floating point sine and 
cosine value of its single precision Hheot ing oint radian argument. 
The call is standard call-by-reference. MTHSGSINCOS_R7 is a special 
routine which is the same as MTH$GSINCOS, except a faster non- 
3 Stendere 358 call is used with the argument in RO and no registers 

; are saved. 


MTHSGSIND and MTH$GCOSD are functions which return the weet ing point 
sine or cosine value of their single precision floating point argu- | 
; ment (de regs). The call is standard call-by-reference. 
: oye « R? and MTH$GCOSD_R7 are special routines which are the same 

; aS MTHSGSIND and MTHSGCOSD” except a faster non-standard JSB call is 

; used with the argument in RO and no registers are saved. 


; MTHSGSINCOSD is a routine which returns the ese point sine and 
; cosine value of its single precision floatin fers egree argument, 
; The call is standard call-by-reference. MTHSGSINCOSD_R/ is a special 
; routine which is the same as MTH$SGSINCOSD, except a faster non- 
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39 | ; stenderd J88 call is used with the argument in RO and no registers 
; are saved 
6 § 
Z 3; VERSION: 1 
¢ t HISTORY: 
66 ; AUTHOR: 
° : MARY PAYNE & JUD LEONARD, 25-MAY-78: Version 0 
$3 ; MODIFIED BY: 
9 4 1-1 Tryggve Fossum, 28-May-78 
3 HF wenn 
0 75: 
4 16 3; HISTORY: 
oo | TON ng nanek, 25-0UNo78: version 2 
5 » 25-MAY-78: n 
$00 8 3 nines 
44 ? : MODIFIED BY: 
99 g : 2-004 Jeffrey C. Wiener 1-APR-83 
4 : : Edit history for Version 2: 
00 86 : 2-001 - Orig 
000 : 2-002 - change nO THSAL 4_OV_PI to MTHSAL_4_OV_PI_V. RNH Oe tet 
444 a8 : -003 - Rod Ped REDUCE” LARGE to correct Bug reported in QAR 
4 90 :. 2-004 - This. code is a modified version of MTH$DSINCOS. The need to 
00 91 ; modify the previous version of MTH$GSINCOS was brought about by 
0000 35 3 the needs of microVAX. G_floating routines are_not permitted to 
0000 93 ; use H_floating point instructions. JCW 1-APR-83. 
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VX 


X$GSINCOS : § 
petty EC 


0000C16C 


00000000 0000010 
20185444 21FB4009 
B2BBBECC 463A4035C 
21D27F33 D97C4022 
385E2955 6A7A402F 
A775E9BB FDBB4035 
C8826DC9 SF 304004 


00000000 80004066 
00000000 8000C066 
C1F81A63 ASDC3EBC 
6C1716C1 C16C3FB6 
D3912529 F46A3F 7D 


. U 
G_2_0V_PI: 
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-SBTTL DECLARATIONS - 


INCLUDE FILES: 
EXTERNAL SYMBOLS: 
-DSABL GBL 
«EXTRN MTHSAL_4 OV_PI_V 
.EXTRN MTHSSSTGRAL 


-EXTRN MTHSK_FLOUNDMAT 
-EXTRN MTHSSJACKET_TST 


EQUATED SYMBOLS: 


14:28:86 


MTHSJACKET.MAR 


X_1_0V_45 = *xC16C 

MACROS: 
SSFDEF ; Define SFS (stack frame) symbols 
SPSLDEF ; Define PSL$ symbols 


PSECT DECLARATIONS: 
.PSECT _MTHSCODE 


; OWN STORAGE: none 


CONSTANTS: 


G_M1: 


QUAD %*x000000000000C010 


G_P1_OV_4: 
-QUAD *X2018544421FB4009 
G_9_PI_OvV_4: 
QUA “XB2BBBECC463A403C 
G_3_P1_Ov_4: 
.QUA *X21D27F 3309704022 
G_5_P1_0v_4: 
. GUA *%385E29556A7A402F 
G_7_P1_OV_4: 
“XA775E 9BBF DBB4035 


QUAD = “*XC8826DC95F 304004 


45: 
nes -QUAD %*x0000000080004066 
“x000000008000C066 
G_SMALLD: 
QUAD “XC1IFBIA63A5DC3EBC 
6_1_0v_45: 
sQUAD *X6C1716C1C16C3FB6 
G_CONVERT: 


“XD3912529F 46A3F 7D 


PIC,SHR,LONG,EXE,NOWRT 
3 program section for math routines 


; =1 

; 0.7853981633974483E+00 
; 0.7068583470577035E+01 
; 0.2356194490192345E+01 
; 0.3926990816987241E+01 
; 0.5497787143782138E+01 
; 0.6366197723675813E+00 


; 0.4500000000000000E +02 
; ~.4500000000000000E +02 


0.4268868231257969E-06 
0.2222222222222222E-01 


; 0.1828292519943295E-02 
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wu bOe BECLARATIONS - peclatetins Part of Modu 6-S =}382 4: 3 6 LCMTHRTL.S REI UVXGSINCO.MAR; 1 ° (2). 
1 § G_90_OV_PI: 
C1FB81A63 ASDC405C ! i « “geasil ites *XC1FBIA63A5DC405C 3 0.2864788975654116E +02 
C1FB81A63 ASDCOO6C $ : 5 -Q0AD “*XC1FBIA63A5DC006C 3 0.3187183529933799-306 
he 
7 1 4+ PI_OV_2 
07 1 pi/2 
0000500 eiresl? 07 160 -QUAD *x000050 oe reno}? 3 0.1570796326794897E+01 
4 Os E71 t 161 » QUAD mt ye rtiae | 
0 A6265(B1 166 ~QUAD “*X A § ¢B1 : 0.5721188726109832E-17 
A 44 ye 16 QUAD aA8 8 A F Bye 
WO0A 7544 ‘ 164 QUAD “Xx A0007544396C ; 0.4335905065061890E-34 
008 24E0357C 4 192 ‘ -QUAD %*x0000880024E037C0 
£0002 008 #)resd 9 OOA 169 as QUAD Risse et Bg htt : 0.3141592653589793E+01 | 
oie 90% E81 OOA 168 QUAD ms Hh Ec E81 
00003000 g 8 C8A 008 192 QUAD %*x000030 $ 33C8A : 0.1144237745221966E-16 
QOOOEO0O SCO63AC4 008 170 QUAD “*XOOQ00EQ00SCO63AC4 
aig 94 As art 00C 171 QUAD male's 44 938FC : 0.8671810130123781E=-34 
00005300 27043744 id \e QUAD %x0000530027043744 
0D0 174 ; 3epi/2 
00007800 Der ca0se 000 86175 UAD 2400007800097 4032 3: 0.4712388980384690E+01 
00004800 ccs 393 008 176 QUAD “xX ot SHES 3E9C 
00004800 79393CCA O0EO 177 QUAD %*X0000480079393CCA : 0.3834758505292833E-16 
00002800 (Ace rece Boe e 178 QUAD Mets se tert es +44 
sig 00 Z3c bene OF 179 QUAD m4 0059(D396 ; 0.1300771519518567E-33 
00005000 6£A137A0 bore 130 heal QUAD *x000050006EA137A0 
00005000 21FB4039 0100 136 , QUAD %*x0000500021FB4039 : 0.6283185307179586E+01 | 
100 10B43E91 010 18 QUAD mee Sh yea H | 
0000 si + C9A 011 133 QUAD %*x000030006 SECA : 0.2288475490443933E-16 
0000E 63AD4 011 185 QUAD a AOR DE DDR ToS AD4 
80000 $00 0129 705 120 186 QUAD mets 4 zee : 0.1734362026024756E-33 
00005300 704437 0128 187 QUAD %x0000530070443752 
0130 3=—s:188 | 
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XSGSINCOS Floating Point Sine, Cosine od alice 1 SEP=1984 05:25 VAX/VMS Macro v04-00 Page 
pati ATHSESINCOS = Radian arguments mist 7 99:93:82 UMTMRTL. SREJUVXGSINCO.MAR: 1 . 4 
F -SBTTL MTHSGSINCOS - Radian arguments 
7 > FUNCTIONAL DESCRIPTION: 
; i§ : The GSIN, GCOS and GSINCOS routines are based on octant reduction. Given an 
é Y : argument. x, it is written in the form 
5 : : x = 118(2epi) + Ie(pi/4) + ¥1, 
3 4 : where 11 and I are intggerss 0 =< I < 8 and 0 =< Y1 < pi/4. Since GSIN and 
> 5 ; GCOS have a period of 2*pi it follows that 
5 : : GSIN(x) = GSIN(I®(pi/4) + Y1)) _and 
2 3 : GCOS(x) = GCOS(I*(pi/4) + Y1)) 
035 0; Using the trigonometric identities for the sum and difference of two angles, 
8 38 1 ; the following table can be generated: 
0 3 5 : If I = then GSIN(x) = and GCOS(x) = 
0338 5: 0 GSIN(Y1) G6COS(Y1) 
0358 8 3 1 GCOS (pi/4-Y1) sriieicevt 
0358 ; : ¢ GCOS(Y1 -GSIN(Y1) 
Ba28 a3 stints -Y1) sayasiere* -Y1) 
0358 39 ; 4 -GSIN(Y1) -GCOS(Y1 
0358 40 ; 5 Sar -GSIN(pi/4-Y1) 
0358 41; § -GCOS(Y1) GSINC(Y1) 
baep 4g : -GSIN(pi/4-Y1) GCOS (pi/4-Y1) 
0358 44 5 Let Y be defined as Y = Y1 if I is even and Y = pi/4 - Y1, if 1 is odd, then 
0358 45 ; each entry of the above table is of the for +/-GSIN(Y) or *+/=GC0S(Y). Based 
b32 +8 ; on the above remarks, the GSIN, GCOS and GSINCOS routines process the input 
035 47 ; argument x, to obtain I and Y, and based on | selects a suitable polynomial 
B32 rt : approximation, p(y), to evaluate the desired fuction. 
035 50 ; 
Bee 3) ; INPUT PARAMETERS: 
00000004 b32 3§ . LONG =4 
00000004 035 54 x = 1*LONG ; x is input angle in radians 
00000008 035 55 sine = ctor : sine is GSIN(x 
0000000C 8 28 2$ cosine = 5*LONG 3; cosine is ats tx) 


ueeppe ines 
-004 


; Return cosine of argument 


OOF C -ENTRY MTHSGCOS, 


MTHSFLAG_JACKET 


“M<R2, R3, R4, RS, Rb, R7> 


PLES PWN 
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Floating Point Sine, Cosine and Sincos 1g-SE 5 VAX/VMS Macro V04-00 Page Vv) 
ATHSGSINCOS = Radian arguments 6-5 Pa} 984 4: 93; 3 MTHRTL.SRCJUVXGSINCO.MAR; 1 . | =( 
59 | 
$4 
$3 ; Return sine and cosine of argument 
: 
OOFC . 96 -ENTRY MTHSGSINCOS, “M<R2, R3, R4, RS, R6,R7> 
3A 68 MTHSFLAG_JACKET 
6D 00000000 ' GF 9E 5A MOVAB G*MTHSSJACKET_HND, (FP) 
61 3; set handler address to jacket 
61 3; handler 
61 
61 $8 
50 4 BC 50FD 61 0 MOVG ax (AP) 
09go3e6 ‘fF 16 0366 371 JSB atnsesiNncos R? 
08 D $95 i mMOva RO, a@sine (AP) 
0c BE 52 7D 7 MOVQ R2, @cosine(AP) 
04 0374 74 RET 
0375 75 
0375 76 
0375 77 
0375 78 -SBTTL MTHSGSIN 
bare 79 
0375 80 ; 
0375 81 ; Return sine of argument 
b372 Ht 3 
0375 
bare 4 
OOF C 3 3 $2 -ENTRY MTHSGSIN, “M<R2, R3, R4, RS, Rb, R7> 
0377 87 MTHSFLAG_JACKET 
6D 00000000'GF 9E 0377 MOVAB G*MTHSSJACKET_HND, (FP) 
037E ; set handler address to jacket 
037E : handler 
f 
50 04 BC SOFD M3 MOVG @x (AP) 
O00004FB'EF 16 038 JSB MINSGSING , 
04 bh RET 
A 
. 
; . -SBTTL MTHS$GCOS 
8 A 
A 
A 
A 
A 
A 
C 
C 
C 
| 
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00000000 ' GF 


08000584" EF Hg ast 


| 


MOVAB 


MOVG 
JSB 
RET 


G*MTHSSJACKET_HND, (FP) 


: set handler address to jacket 
3; handler 
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gecei" Sine, Cosine ond elnews a abs et Og: 93: §? YA/Nes Be v04 


00 P 1 
D = Degrees 6-SEP-1984 on ( 


MTHRTL. SREJUVXGSINCO.MARS 1 


nw 


tin | 
INC | 
9F 411 -SBTTL MTHSGSINCOSD - Degrees 
OF tig 
oF 41 
SF 414; 
44 rE ; FUNCTIONAL DESCRIPTION: 
oF 419 : The GSIND, GCOSD and GSINCOSD routines are based on octant reduction. Given | 
ar $13 3; an argument, x, it is written in the form 
oF 4 0 : x = 11#360 + 1845 + Y1, 
OF 64 § : where 11 and | are integers =< 1 < 8 and 0 =< Y1 < 45. Since GSIND and 
oF 4 33 GCOSD have a period of 360 te” follows that 
9F 405 ; GSIND(x) = GSIND(1*45 + Y1) and 
$e ? $ 3 GCOSD(x) = GCOSD(1*45 + Y1). 
O39F 428: Using the trigonometric identities for the sum and difference of two angles, 
baoF ; 4 3; the following table can be generated: 
O39F 431 : If l= then GSIND(x) = and GCOSD(x) = 
039F 235 ; 0 GSIND(Y1) GCOSD(Y1) 
O39F 434; 1 Ba bse Sea) GSIND(45-Y1) 
Baar 435 ; § GCOSD(Y1) -GSIND(Y1) 
O39F 436 ; GSIND(45-Y1) -GCOSD (45-1) 
O39F 437 ; 4 -GSIND(Y1) -GCOSD(Y1) 
O39F 438 ; 5 -GCOSD(45-Y1) -GSIND(45-Y1) 
039F 439 ; 6 -GCOSD(Y1) GSIND(Y1) 
Baer rr ; 7 ~GSIND(45-Y1) 6C0S(45-Y1) 
039F 226 : Let Y be defined as Y = Y1 if I is even and Y = 45 - Y1, if I is odd, then 
O39F 443 ; each entry of the aheve table is of the for +/=GSIN(Y) or +/=GCOS(Y). Based 
O39F 444 ; on the above remarks, the GSIND, GCOSD and GSINCOSD routines process the input 
O39F 445 ; argument x, to obtain I and Y, and based on I selects a suitable polynomial 
Baer 268 ; approximation, p(Y), to evaluate the desired fuction. 
039 res 
00000004 9F 449 LONG = 4 
00000008 te 450 sind = pig ee: 
0000000C O39F 451 cosd = 3*LONG 
O39F 452 


M11 
VXSGSINCOS Floating Point Sine, Cosine and Sincos 16-SEP-1984 02:05:25 VAX/VMS Macro v04-00 Page 1 
wu bOe ATHSGSINCOSD - Degrees aay 1 96:93:82 MTHRTL.SRCJUVXGSINCO.MAR; 1 ° 18) 
OOFC +4 $38 eENTRY MTHSGSINCOSD “M<R2, R3, R4, RS, Rb, R7> 
Al 456 MTHSFLAG_JACKET 
6D Q0000000'GF 9E Al MOVAB G*MTHSSJACKET_HND, (FP) 
A 3 set handler address to jacket 
3; handler 
A 457 
50 04 BC SOFD A 458 MOVG aXx(A 
OOOO061F'EF 16 O3AD 459 JSB Rywbgs{ncosp, R7 
Be BC 9 7D B $69 MOVQ ay. & 
C BC 7D B 461 MOVQ ee ba 
8 BB $06 
04 BB 46 RET 
03BC 464 
038C 465 
038C 298 
OOFC tts 46 ~ENTRY MTHSGSIND “M<R2, R3, R4, R5, Rb, R7> 
O3BE 468 
O5BE 469 MTHSFLAG_JACKET 
6D 00000000 ' GF 9E tra) MOVAB G*MTHSSJACKET_HND, (FP) 
03C¢ 3; set handler address to jacket 
beta 3; handler 
03C5 
03¢5 470 
50 04 BC SOFD O35 471 MOVG ax (AP) 
0000067D'EF 16 O3CA 47g JSB MIHSGSIND. R? 
0300 47 
04 0300 474 RET 
03D1 475 
03D1 476 
0301 477 
OOFC 03D1 $78 -ENTRY MTHS$GCOSD “M<R2, R3, R4, RS, Rb, R7> 
0303 «47 
§303 480 MTHSFLAG_JACKET 
6D 00000000 ' GF 9E 0303 MOVAB G*MTHSSJACKET_HND, (FP) 
O3DA 3 set handler address to jacket 
O3DA : handler 
O3DA 
O3DA = 481 
Rapa $86 MOVG @X(AP), RO 
DF 4 JSB MTHSGCOSD_R7 | 
Beee ret 
3€5 = 485 
03E6 486 
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ting Poi 
INCOS_R7 P-1984 8:56 CMTHRTL.SRCJUVXGSINCO.MAR; 1 (10) 
ES iB -SBTTL MTHSGSINCOS_R7 
E 490 ; This routine computes the GSIN and GCOS of the G-format value of RO/R1. The | 
E6 491 ; computation is performed one of three ways depending on the size of the 
Ee 138 : input argument, X | 
6 494: 1) If txt < pi/4, then X is used directly in polynomial approximation 
—E6 495; of GSIN and GCOS. 
E6 £36 3 2) If pi/4 =< ixi < 9epi/4, then the subroutine REDUCE_MEDIUM is called 
6 497 ; to reduce the orguaent to an equivalent eroynens in radians, Y, ane 
E6 $38 3 the octant, I, conta anine the Srqusents is then evaluated in 
ES 499 : polynonials ¢ chosen, as a unction : on compute GSIN(X) and GCOS(X). 
E 280 : 3) =< iXi, then the ereuttes REDUCE_LARGE is called to 
E6 501 ; poduce tye argument to an equivalent Bein in cycles, Y ane the 
O3E6 206 ; octant, contanining the argument. Y is then evaluated 
baee 207 3 nd pd chosen as a function of I, to compute GSIN(X) ne GCOS(X). 
b3e8 505 MTHSGSINCOS_R7:: 
© 30 NR Osea 809 Bees og SgNCOS ‘pede | 
000003F7'EF 16 O3EC 508 JSB SINTOS : RO/R1 = GSINCiX!), R2/R3 = GCOS(X) 
50 50 52FD O3F2 509 MNEGG RO, RO ; RO/R1 = GSIN(X) 
05 O3F6 510 RSB 
te 31) SINCOS: 
56 8000 8F AA O37 318 ae since #*xX8000, R6 : R6/R7 = IX! 
56 FCO7 CF SIFD O3FC 515 ~  CMPG  G.PI_oV ; Compare pi/4 with {X! | 
35 14 0402 516 BGTR Smale Sthéos: ; No argument reduction is necessary 
56 FCO7 cr 51FB et af CArG -PT_OV_4, R6 3; Compare 9*pi/4 with {Xx! 
ooc2 31 O40C 519 BRW heat, ‘sianes ; Use special logic for IX! > 9*pi/4 
Odor 831 ; 
O40F $33 > pi/4 =< IX! < 9epi/4 
00000720'EF 16 0408 324 is: JSB REDUCE_MEDIUM : Medium argument reduction rout ine 
3 = reduced argumen 
0415 556 : R2 = octant . 
. 2. See bei2 $37 MOVQ R4, -(SP) ; Save reduced argument on stack 
7E 56 7D 418 528 MOVQ R6, -(SP) : 
52 DD 0418 44 PUSHL R2 3 Save, ectant bits on stack 
000005A8'EF 16 pele 530 JSB M_COS ; RO/R GCOS (Xx 
se 0 Be 531 MOVL (SP)+, R2 ; R2 = ‘mastet bits 
56 E D 0426 236 MOVQ (SP)+, R6 : 
54 ef 7D 0429 23 MOVQ (SP)+, R4 ; R4/R7 = reduced argument 
7E 7D 042C 534 mOVa RO, (SP) > Save GCOS(X) on stack 
Q0000526"EF 16 O4eF 535 JSBMLSIN : RO/R1 = GSIN(X) | 
52s 7D 0455 238 MOVQ (SP)+, R2 : R2/R3 = GCOS(X) | 
05 " 3 53 ‘ RSB 
239 335 : Logic for small arguments. ixXi < pi/4. 
ye 
56 4000 8F BI £33 ste ates Sari #*xX4000, R6 Compare 1/2 with ‘Xi 
. e ia ' 
5 15 3e 544 BLEQ 2$ : Sufticent overghang not available 
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FCES CF 9 8 
4 6 

55 _FFFFO7FF 8F 
7 ¢ rr 
5 6 6 
52 6 

0 
52.10 

50 ef 

Sh 4 

54 =10 

54 9? 

50 4 

52 oF 

6— 5 

FCEB CF 38 32 
0 56 

50 0-56 
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nt Sine, Cosine and Sincos 16-SEP-1 


-SEP=1986 11:58:86 


#*X3E60, R6 

Is 

R + Gg, RG 

Ro. ACOSLENRI- ¥ COSTAR! 
(SP) 

R isp sy 

RG, #SINLENR=1, SINTBR 
R6, R 


Da Be Be Be Be Be Be Se Se Se Se 
~ 
@D 


R6, R 
(SP)+, R2 


R2/ 
#1.0, R2 : R2/R3 = 1.0 = GCOS(X) 
#18600" RO + RO/RI = 1K! 
R6, R6, RS R4/RS = .x*2 


AX/VMS Macro V04-00 P 
MTHRTL.SRCJUVXGSINCO.MAR; 1 


Compare with 2*-27 
ae Te ae evaluation is needed 


X*X 


R4, -(SP) 
R4, #COSLENR2-1, COSTBR2; ae = pei 


- + RG 

#°XFFFFOTFF, RS : 

R4, R6, -(SP) : ($P) = XL LO. Last 2 

R6, RO. R2 + R2/R3 = X + XH 

(SP)+, R : 4 /R3 = XLO®(X + XHI) 
$ ; Check for A r, = 0 
#°x0010, R2 + R2/R3 = A2/ 

R2, RO : RO/R1 = Q(Y*2) = A2/2 
RG. RG : RG/RS = KHI* 
#*%0010, R4 + R4G/R5 = XHI*2/ 
#1, R46 : R4/R5 = XHI*2/2 = 1 
RG RO : RO/R1 = GCOS(X) 

(SP), R2 t R2/R3 = X*2 
RO, (SP) : Save GCOS(X) 
R2. #SINLENR=-1, SINTBR : RO/R1 = Q(xX*2) 
R6, RO + RO/R1 = X#Q(X*2) 
R6, RO : RO/R1 = GSIN(X) 
(SP)+, R2 > R2/R3 = GCOS(X) 


REDUCE _LARGE ; RG/R7 = 


; RO/RI = 


. 
zz 
= 
Ree Oe Oe Oe Oe Be Oe Oe Oe ee 


R2 = 
R2/R3 = 
RO/R1 
R2/R3 


"2G GCOS (X) 


: Reduced 6 orqugent 
octant bits 


GCOS 
GSIN 
GCOS 


(x) 
(x) 
OS(X) 


pave GCOS(X) on stack 


RO/R1 = + ae 

RO/R1 = rt 

RO/R1 = GSIN(X) 
R3 = GCOS(X) 


: /R 
; R4/RS = XHI, big Me A ble bits 


reduced argument (in cycles) 
; R2 = octant bits 

; Save octant bits on stack 

; weve. reduce 
ent on stack 
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50 53FD 

10 «18 

seen Yh 16 
50 0 52FD 

05 

50 8000 8F AA 
50 =FAF3 CF aire 
50 =-FAFS gf 51FD 
8 19 


00000720°EF 16 


vr" of BS 


50 3660 of °) 
56 5 70 
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.SBTTL MTHSGSIN_R7 


ro V04-00 
MTHRTL.SRCJUVXGSINCO.MAR; 1 


4 ; This routine computes the GSIN of the G-format value of RO/R1. The 


; computation is performed one of three ways depending on the size of the 


input 


Sete Ge Ge Se Se Ge Se Ge Se Se Se 


MTHSGSIN 


SIN: 
POS_SIN: 


argument, Xx 


Page 


1) if HY < pi/4, then X is used directly in a polynomial approximation 


) N. 
2) If pi/4 =< txt < 9epi/4, then the subroutine REDUCE_MEDIUM is called 


to reduce the argument to an equivalent argument in radians, Y, a 
the octant, I, containing the argument. Y is then evaluated in a 
pg znoe el chosen as a function of I to compute GSIN(X). 

3) If Sepi/4 =< IX!, then the subroutine REDUCE _LARGE is called to 


reduce the argument to an equivalent argument in cycles, Y 
octant, I, contanining the argument. Y is then evaluated in a 


and t 


polynomial chosen as a function of I to compute GSIN(X). 


R7:: 

TST1G RO 
BGEQ POS_SIN 
JSB N 
MNEGG RO, RO 
RSB 


BICW2 #*x8000, RO 
G_PI_Ov_4, RO 

BGTR SMALC_SIN 

CMPG G_9_PT_Ov_4, RO 

BLSS LARGE_SIN 

=< IK! < 9epi/4 


JSB REDUCE MEDIUM 


CASEB R2, #1, #7 


“WORD P7COS"R-1$ 
“WORD N7SINTR-1$ 
“WORD N7SINTR=1$ 
“WORD N7COS"R-1$ 
“WORD NTCOS7R-1$ 
“WORD P7SINTR-1$ 
“WORD © PSIN7R-1$ 


: Logic for small arguments. iXi < pi/4. 


SMALL _SI 


N: 

CMP #*X3E60, RO 
BGEQ 1$ 

MOVQ RO, R6 


Check the sign of RO 


RO/R1 = GSINCiX!) 
RO/R1 = GSIN(X) 


RO/R1 = IX 


Compare pi/4 with iX! 

No argument reduction is necessary 
Compare 9*pi/4 with IX! 

Use special logic for iX! > 9*pi/4 


Medium argument reduction routine 
R4/R7 = Y = reduced argument 

R2 = octant 

Branch to one of four polynomial 
evaluations depending on the 


octant bits. 


; Compare with 2*-27 


No pel ynentat evaluation is needed 
R6/R7 = X 


nd 


he 
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; Floati 
HS$GSIN_ -SEP=1984 MTHRTL.SRCJUVXGSINCO.MAR; 1 (11) | 
0 44FD 44 659 MULG2 nt ; . /R1 = 
FC61 CF 6 7728 re 268 PUL Gg R Bane SINTBR ; ss Yd : qix"2)_. 
; x 
6 sore 5 ; 066 ADDG Re. ; RO/R1 = GSIN(X) 
0 : 663 1$: RSB 
+e 
3 096 LARGE_SIN: 
QOOOO7C2"EF 16 22 pots JSB REDUCE _LARGE 3 erat = reduced argument (in cycles) 
$ = octan s 
54D 55 96? L_SIN: TSTL R4 ; Check for degenerate case 
14 1 26 6/9 BEQL DEGENERATE_CASE_SIN 
07 OO S52 &F 2° 6 CASEB R2, #0, #7 
0771' 0566 674 1$: «WORD P_SIN_C-1$ 
0689" 056 675 «WORD P_COS_C-1$ 
p68? 56A 676 «WORD P_COS_C-1$ 
1° 056C 67 «WORD P_SIN_C-1$ 
0767" 056 678 «WORD N_SIN_C-1$ 
0708" $70 679 «WORD N_COS_C-1$ 
14, B b25 680 «WORD N_COS_C-1$ 
0767* 0574 681 «WORD N_SIN_C-1$ 
Bit 
B25e one DEGENERATE_CASE_SIN 
52 8A 0576 686 BICB2 #1, R2 ; Compute index as (R2 = 1)/2 
52 52 FF 8F 9C 0579 687 ROTL #1, R2, R2 
3.©6—SfC«O0 8F 0386 tt) CASEB- R2, “#0, "#8 
0868" $38 690 1$: «WORD P_ONE-1$ 
087D' bege 691 -WORD UNFL -1$ 
0870' B2 $ O36 «WORD N_ONE-1$ 
087D' beee 69 ~ WORD URFL -1$ 


o 
oOo 
=> 
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XSGSINCOS : Flo 7 
petty rude db) 


-SBTTL MTHSGCOS_R7 
This routine computes the GCOS of the G-format value of RO/R1. The 


; computation is performed one of three ways depending on the size of the 
; input argument, X. The processing is the same as described for MTHS$GSIN_R4. 
MTHSGCOS_R7:: 
30 53FD TSTG RO : Check for reserved operand 
2° P0090 F AA BicWw2 #*x8000, RO 3 RO/R1 = IX! 
GQ FA71 CF SIFD CMPG G_PI_ OV_4, RO : Compare pi/4 with {X! 
2. 14 BGTR SMALC_COS ; No argument reduction is necessary 
50 FA71 CF S1FD CMPG G_9 PT_Ov_4, RO 3 Compare 9*pi/4 with {x! 
48 «(19 BLSS LARGE _Cos ; Use special logic for ix! > 9*pi/4 


: piss =< 1K! < 9epis4 


Q0000720°EF 16 JSB REDUCE _MEDIUM 


Medium argument reduction routine 

R4/R7 = Y = reduced argument 

R2 = octant 

Branch to one of four polynomial 
evaluations depending on the 


~~ RR RF M_COS: CASEB R2, #1, #7 


0648' 1$: .WORD N_SIN_R-1$ 

0648" “WORD N7SINTR-1$ 

OSEF' [WORD N7COS"R-1$ 

OSE “WORD N7COS7R=-1$ 

06 4 “WORD PTSINTR-1$ 

0652" “WORD PTSINTR-1$ 

059C" “WORD P7COS"R-1$ 
059C" “WORD P7COS7R-1$ octant bits. 


; Logic for small arguments. iX! < pi/4. 


SMALL_COS: 
CMPW 


ewan eal eal eal eal ealeal eal ealvalealvalvalvabrabeabrabr ap abrabrabrabvalvabvalvabvalvalvalvalvalvalvalvalvalvalvalvalval val al al alval aval Av al Al ah abv ab al aM mde 

MMMMMMMMMoOCOCVTIOOCOCOCWWWVIVIVIOI@OIVIOIOIOWr\r rr rrr rr Kr PrP OOOWWDWDOOWDODOOWwom Z-—- 

SVT FT MWOVNOM PF VMVS AW NOOOOOO SWS LMM OOWWOWrynyyunwvorForvwery?>y YY YY Yr YrYrYr>? i 
ENO OOO 
FEES FWWIWWIWIWIWIWIWIUIRININININININNNININ 2 2 3 OH HF MHQOODOOOOCOOOCOOOOO 
FPWR OD ODNA UE WN 9 OONO UE WIN $0 OO NAME WWIN  O OONA UE WN O ODN 


50 4000 8F_ BI #*x4000, RO ; Compare 1/2 with ix! 
1 14 BGTR 1$ : Sufficent overghang is available 
-— 6 7D MOVQ RO, R6 3 R6/R7 = : 
57 _FFFFO7FF 8F CA BICL @*XFFFRFO7FF, R7 ; R6/R7 = XHI, 26 most signif bits 
54 50 6 43FD SUBG R6, RO, R 3 R4/R5 = XLO, 27 least signif bits 
95 0 631 BRW NEEDS bouBL E ; Use special logic to obtain overhang 
50 3660 8F 861 1$: CMPW #*X3E60, R ; Compare with 2*-27 
. 2 on BGEQ 3; No priprearet evaluation is needed 
50 0 44FD MULG2 RO,RO ; RO/R1 = X* 
FB49 CF «607 )=— 550 55FD POLYG RO, #COSLENR1-1, COSTBR1; RO/R1 = GCOS(X) 
05 Sk RSB 
50 08 SOFD 747 2%: MOVG #1.0, RO ; RO/R1 = 1.0 = GCOS(X) 
05 74 RSB 
74 
£39 
751 LARGE_COS: : 
OOOOO7C2"EF 16 752 J REDUCE _LARGE ; R4/R7 = reduced argument (in cycles) 
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XSGSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 5: AX/VMS Macro V04-00 Page 19 | 
peti ATHSGSINCOSD R7 a ety 4: 93 3 MTHRTL.SRCJUVXGSINCO.MAR; 1 . (14), 
sit 4 9 -SBTTL MTHSGSINCOSD_R7 
Ps § ; This routine computes the GSIND and GCOSD of the G-format value of RO/R1. 
61F 785 ; The computation is performed one of two ways depending on the size of the 
pif if : ; input argument, X 
6i-F 7 § s 1) If txt < 45, then X is used directly in polynomial approximation 
61F 787; of GSIND and GCOSD. 
61F 788 ; 2) If 45 =< txt, then the subroutine REDUCE_DEGREES is called fe reduce 
61F 789; the argument to an equivalent argument in degrees, Y, and t 
61F 790; octant, I, containing the argument. Y is then evaluated A. twe | 
61F 791; pot ynosie s chosen as a function of I, to compute GSIND(X) and 
Se 3g | 
OIF 794 MTHSGSINCOSD_R7:: 
ag 7 ee 
50 8000 8F AA 0624 797 BiCw2 #*x8000, RO 3 RO/R1 = {Xi 
00000634'EF 16 be oe JSB SINCOSD : Ad = Se epciatt 
50 50 52FD 062F 800 MNEGG RO, RO 3; RO/R1= -GSIND(iX!) | 
oie oe | 
pose 508 SINCOSD: 
50 FOFF CF SIFD 0634 804 CMPG RO ; Compare 45 to [X! 
4 14 Q63A 805 BGTR Saat SINCOSD ; special processing for small arg 
QOOOOABC'EF 16 B676 aoe JSB REDUCE_DEGREES 3 ey = reduced argument 
: = octan 
— a O63 808 MOVQ R6, (SP) 3; Save reduced arg | 
53 oD Bee 809 PUSHL R53 ; Save octant bits 
Sooner e* er 16 0647 810 JSB EVAL. cOSsD ; RO/R1 = GCOSD(Y) 
5 8E 00 064D 811 MOVL (SP)#, R3 ; R35 = octant bits 
6— 7D 8338 aig MOVQ (SP), R6 ; R6/R7 = requces argument 
6E 50 7D 065 81 MOVQ RO, (SP) ; Save GCOSD(Y) 
OOOO06A0'EF 16 3e2e 814 JSB EVAL_SIND ; RO/R1 = OS INDY) 
52 soBE 7D 065C 815 MOVQ (SP)F, R2 ; R2/R3 = GCOSD(Y) 
05 O65F 816 RSB 
oeeo BB 
660 819 SMALL_SINCOSD: 
, ee 94 820 SUBL2 16, SP $ a locete 4 _Lenguerds on stack 
50 7D 066 821 Ova RO, (SP) : argum 
EF 16 0666 8 : JSB SMALL cOSD ; ROvRI GCOSD( ix!) 
50 7D 06 § MOVQ 0, 8tSP ; Save GCOSD(:x Xi) 
BE 7D 06 4 MOVa = (SP) +, : RO/R1 = argument 
EF 16 bor H 5 JSB SMALL_SIND ; RO/R1 = GSIND(X) 
8E 7D 0679 8 MOVQ (SP)+, R2 ; R2/R3 = GCOSD('X!) 
05 067¢ 8 R 
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This rout ine computes the GSIND of the G-format value of RO/R1. 


computat 
argument, X 


1) If iy 
2) If 4 
octant, 


BGEQ 
JSB 
a 
RSB 
G_SIND: 


The 
on is performed one of two ways depending on the size of the input 


Ry 45, then X is used directly in polynomial approximation 


=< ixi, then the subroutine REDUCE_DEGREES is ag ig" fe reduce 
the Srqugens to an equivalent orgqunent 
containing the ar 
) TR chosen as a func 


MTHSGSIND_R7:: 
TSTG 


#*x8000, RO 


45, RO 
SMALL, SIND 
REDUCE _DEGREES 


R3, #0, #7 
P_SIN_b-1$ 
P"COS"D-1$ 
P~COS"D-1$ 
P"SIN7D-1$ 
N7SIN7D-1$ 
N7COS"D-1$ 
N7COS"D-1$ 
N7SIN7D-1$ 
G_SMALLD, RO 
1$ 

Boe Re _D 

RO 


ty 
25° SMALLEST_DEG, RO 


aE oh R2 


#*X006 
R2 


rage a 


in degrees, Y, an 


is then evaluated - two 
ion of 1" to compute GSIND(X). 
RO/R1 = X 
RO/R1 = GSIND(IX!) 
RO/R1 = -GSIND(iX!) 
RO/R1 = {Xt 


Compare 45 to iXi 

special processing for small arg 
RO/R7 = reduced argument 

R3 = octant 


Compare 180/pi*2*-27 with ‘x: 
No polynomial evaluation is 
necessary 


Check for zero 
Return if RO = 
Check for possible underflow on 
conversion to radians 
Under flow will occur on conversion 
R2/R3 = (pi/180 - 2°-6)*ix! 
RO/R1 = iXi#2*-6 
= GSIND(!X:! )= (pi/180) ix: 
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oint Sine, Cosine and Sincos 16-SEP-19 


4 -SBTTL MTHSGCOSD_R7 


§ ; This routine computes the GCOSD of 
; is performed one of two wa 
; Details are given in the d 


MTHSGCOSD_R7:: 
9 TG 


TS 

9 BICw2 

9 CMPG 

94 BGTR 

95 JSB 

i 

\) 

98 EVAL_COSD: 

99 CASEB 

00 1$: WORD 
901 WORD 
90 ~ WORD 
90 . WORD 
904 ~ WORD 
905 WORD 
90 WORD 
90 ~ WORD 
908 
909 
910 SMALL_COSD: 
911 CMPG 
aig BGTR 
91 MOVaQ 
914 BRW 
915 1$ MOVG 
916 RSB 
917 


RO 
#*x8000, RO 
SRA 


ww 


vwuUszzzzuv0y 
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G_SMALLD, RO 
1$ 


B tos’ p 
#T, RO 


he G-format value of RO. 
depending on the size of the input argumert, X 
ussion on MTH$GCOSD_R4. 


it_coso 
REDUCE_DEGREES 
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Check for reserved operand 
RO/R1 = IX! 
Compare 45 to {Xi 


R6/R7 = reduced argument 
R35 = octant 


Compare 180/pit2*-27 with ix: | 

Check if polyinomial evaluation is 
necessary. 

POLY needed 

RO = 1. = GCOSD(:X!) 


The computation 
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XSGSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 02:05:25 VAX/VMS Macro v04-00 Page 
petit MTH GCOSDR? ae 94:93:83 EMTHRTL. SREIOVXGSINCO.MAR: 1 . (48) 
135 338 
7 921 -SBTTL REDUCE_MEDIUM 
15) 958 ; 
7 924 ; This routine assumes that the absolute value of the argument, X, is in RO/R1 
7 925 ; and that pi/4 =< ix! < 9*pi/4._ It returns a pair of G-format values for the 
7 9 § 3; reduced argument: YHI in R6/R7, and YLO in R4/R5. The octant bits are 
4 4 3 returned in R2. 
7 9 5 ; The reduced organs is obtained by locating the octant that X is in through 
g 4 ? 3 a binary search and then subtracting off a suitable multiple of pi/2 
7 932 ° 
7 6 2 : REDUCE MEDIUM: 
50 8000 8F AA 07 9 934 BICW2 #*x8000, RO 3 RO/R1 = {X! 
50 F8F6 4 site o/ H 332 ys 5_P1_Ov_4, RO : sei ‘os Semtl 
50 FBEG CE SIFD O79R 987 CHPG 3_P1_0V_4, RO re ‘ 
0S 15 07 R 938 BLEQ gs y IX! >= 3epi/4 
ae oe, BF 5 939 MOVL #1, R2 ; First quadrant 
7 «171 07 8 40 BRB SUBTRACT 
52 03 00 O73A ang 3$: MOVL #3, R2 3; Second quadrant 
12. «11 H : 343 BRB SUBTRACT 
50 FBES cr StF O73 945 5$: CMPG G,7-P1_0v_4, RO ; er ee 
Ss ' = 1 
52 05 00 0747 309 MOVL #5, R2 ; Third quearent 
0511 Bree 338 BRB SUBTRACT 
52 07 dO O74C £950 7$: MOVL #7, R2 : Fourth quadrant 
00 11 O74F 951 BRB SUBTRACT 
O78] 988 
O78] 985 S New 
; Ne 
57 5202 C5 0751 986 “  MULL3 #2, R2, R7 
53 FA ATG E 0755 9 MOVAW Ae R3 ; R3 = index into Pl_OV_2 table 
53 F911 CF43_ ODE O75, 958 wey OVAL PI_OV_2CR33, R3 : R3 = pointer into PI_OV_2 table 
DO 0760 960° MOVL. RO, R6 
57 51 FFFFOSEF er cB 768 961 BICL @*XFFFFO7FF, R1, R7 ; R6/R7 = XH = 26 HIGH BITS OF X 
7E 20 43FD 768 96 SUBG3 6, RO, -(SP) : SP = XL = 27 LOW ORDER BITS 
6 3 42FD 07 96 SUBG (RS), R6 
ge : 42FD 0774 964 SUBG (R35)+, (SP) 
6 E 40FD 0778 965 G (SP)+, R6 
oe a9 ; END NEW 
77C «= 9668 gx SUBG3 (R3)+, RO, R6 ; R6/R7 = 1st approximation to YHI 
15 O77C 34 BLEQ 1$ 3 = YHI' 
D6 O77E 970 INCL R2 ; Adjust octant bits 
54 56 8000 BF AB 0780 971 1S: BICW3 #*x8000, R6, R4 ; R4 = high 16 bits of !YHI'! 
54 D10 8F B81 0786 378 CMPW #*x3010, R4 3; Check for at least_6 Signiticent bits, 
? 97 ; ie, CMPW exp of 53-6=47 with (R4) 
7 974 3; so there is at least 6 bits of overhang 
0 14 07 975 BGTR  NOT_ENOUGH_BITS ; 


K 12 
XSGSINCOS Floating Point Sine, Cosine ond Sincos 16-SEP-1984 AX/VMS v04-00 p 3 
hoe REDUCE. MEDIUM COs pee loe $1i5b:8e TATANTy eke dovNccinco.mans1 2! cf 
78D 7 
54 $ 7 07 7 mova Rg, R4 : R4/R5 = YHI' 
D4 07 7 CLRL. OR t R6/R7 = high 21 bits of YHI' = YHI 
4 6 43FD t 7 SUBG R6 Ra : R4/ Rd = low bits of YHI' 
4 634 FD p SUBG2 = (R48), + R4G/R5 = YLO 
> hon 
ao ee ae 198 NOT -ENOUGH BITS: 
79¢ 5; BICL2 #*x FD ; 
50 00070000 cA te : Bi cL eset $08" RO : 
54 83 43FD O7A 8 SUBG RD. 8 (R3)+, RG ; 
é 42FD O7A SUBG + R6/R7 = YHI 
4 é 40FD 0780 990 ADDG Red)” ons + R4/R5 = -YLO 
54 8000 8F aC 7B4 931 XORW  #*x8000, R4 + R4/R5 = YLO 
05 0789 36 RSB 
7BA 99 
7BA 994 
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VX$GSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 02:05:25 VAX/VMS Macro V04-00 Page 24 | N 
peti REDUCE. LARGE g-SEp=19 4 99:03:88 MTHRTL.SRCJUVXGSINCO.MAR; 1 ’ (95) | ~( 
7BA -SBTTL REDUCE_LARGE 
ae 
7BA : This routine is used to reduce large arguments (iXi >= 9*pi/4) modulo pi/4. 
7BA ; It returns the reduced argument, Y, in R4/R7 in units of cycles, and returns 
yea ; the octant bits, I, in R2. 
7BA : The method of reduction is as follows: 
7BA : x*(4/pi) = Zanetsca/pi) here n is an integer and 1/2 =< f < 1 
7BA ; = 2*(n=-53) #( a53e¢)#(4/pi) 
7BA : = (2°53#f)#(2*(n-53) *4/pi) 
7BA ; = K*C, where K = 2°53*f is an integer and C = 2*(n-53)*4/pi 
07BA 3 Let L = K*C modulo 8, where 0 =< L < 8, and let I = the integer(L) and 
O7BA : h = fract(L), then if I is even Y = h, otherwise Y = 1-h 
O7BA : CONSTANTS: 
00000300 07BA L_INT WEIGHT = *x03D0 3; weights exponent by 4 
00000200 O7BA W_TERM_WEIGHT = *x0200 ; weights exponent by 32 
00000400 07BA W-MAX_QEIGHT = *x0400 3 maximum unbiased exponent 
00000386 07BA Ost = *X3B6 ; Used to locate binary point in 
;  MTHSAL_4_OV_PI_V table 
G_2_T0_32: 
00000000 00004210 QUAD *X4210 3; FR 


REDUCE _LARGE : 


The first step_is to_obtain the location of the binary point in the represen- 
tation of C = 2*(n-53)*(4/pi) in two parts - the number of longwords from 

the start and the number of bits from the most significant bit of the next 
lLongword. Also K = 2*53*f must be obtained. 


BICW2 Mh gag RO RO/R1 = iX! . 

ROTL #-4, RO, R3 Shift exponext field 4 bits right 

SUBW2 #w_AbdJusT, R3 Unbias exp and adjust for leading 
zeroes. R3 = location of binary 


SSNS 


50. 8000 8F 
50 FC 8F 
53 0386 8F A2 


oint 
54 53 fe 8F ROTL #-3, R3, R4 Divide R3 by 32 and mull by 4 to get 
54 FFFFFFOS 8F CA BICL2 #*XfFFFFFO3, R4 R4 = # of Longwords (in bytes) fo 
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7 
7 
7 
7 
7 
7 
7D 
70 binary point. 
52 QO000000°EF ODE 43 MOVAL MTHSAL_4_OV_PI_V, R2 tet pese address of MTH$AL_4_OV_PI_V 
able 
52 Soooppeo "Er CO Q7E ADDL WT NSAL 4.0V_PI_V, R2 R2 = address of MTHSAL_4_OV_PI table 
52 4 C2 O7E SUBL R4, R R points to 1st quadword of interest 
53 €0 8F' BA 43 BICB #*XE0, R3 R3(7:0) = # of bits within longword 
50 7FFO 8F AA O7F BICW2 #*X7FFO, RO ; Clear exponent field 
50 4150 8F a8 O7F BISW #*x4150, RO ; RO = 2*21f 
50 4AFD rs CVTGL RO, R 3; RO = High $| bits of K 
51 51 10 9 ROTL #16, R1, RI : R1 = Low 32 bits of K 
9 1 4 BGEQ 1$ ; Check for high bit of R1 set 
D6 : INCL RO ; Adjust RO if R1 is negative 
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xz: 
m 
om 
c-— 
oo 


it 
in 


in computing x*4/ consist of the fir e 
be o 


t 

C and the first 61 fraction bits of C. These bits w 
constant stored in the interger array MTHSAL_4_OV_PI_ 
NOTE: The ASHQ, ADDL, and MULL instructions in the follow sections may 
result in an ergo overflow trap. The overflow incurred is intentional, 
so that the IV bit must be turned off. The IV bit is not restored until 
after all of the necessary fraction bits hav been generated. 


$: MOVPSL =(SP) ; Put current PSL on stack 
BICL #*C<PSLSM_IV>, (SP) ; (SP) = current IV bit 
BICPSW #PSLSM_IV ; Clear integer overflow bit 


mB ee Ge Ge Ge Ge Ge Ge Ge Ge Se 


7E 
6E FFFFFFOF of 


@Woo 
oro 


The necessary calculation to produce the reduced argument can require up to 
nine longwords of temporary work space. This work space will be allocated 
on the stack. The work space will be accessed though the use of three 
registers: R4, R5, and SP. For the purposes of comments the temporary work 
space will be referred to as locations [0 though T8. The stack and its 
pointers will look something Like this: 


<---- SP 


T0: 
TI: 


T2: 
T3: 


T4: <---- R4 


17: 


T8: 


<e--- R5 


The following code allocates the storage and sets up the pointers. 


#36, SP 3; Allocate 9 longuords on the stack 
ADDL #48. SP, R4 ; R4 points to T 
#32, SP, RS ; RS points to T8 


wu 

ws 

we 
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— =) 


; Get C** = €(0):C(1):C(2):C(3) in TS/T8. C(O) though C€(3) are unsigned — 

: integers generated from the binary representation of C. The high three bits 
: ) are the the first three bits to the left of the binary point of C. 

; The remaining bits C(O) and C(1) though C(35) 25 bits to the 

; right of the b nery point of C. Note that the C(i)'s are adjusted to 

; compensate for their signed (rather than unsigned) interpretation in the EMUL 
; instruction. Note also that the representation of C has no more than 15 
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VXSGSINCOS ; Floating Point Sine, Cosine a Sincos Ven 36b=18 5 VAX/VMS Macro V04-00 Page 26 
2006 REDUCE _LARGE SEP=1984 11:28:56 EMTMRTL. SRESOVXGSINCO.MAR; 1 99 4B), 
if 1319 : consequtive ones, so that no carry is possible from the adjustment. | 
1E 1112 ° 
57 54 OC C1 23 1118 ADDL3 #12, R4, R7 ; inital) tae loop counter. R7 points 
67 62 53 79 111? ASHQ R3, (R2), (R7) snife. th the cbfgrer ,queduore so that | 
- ms 4 143 SUBL #4, R7 ; R7 points to 16 
5206 C2 0899 1178 28: SuBL a4, Re : R2 points to next quadword in | 
67 62 53 79 $8 C yt y ASHQ R3, (R2), (R7) ont tt ayadvard” so that a Pe is in 
03 18 $8 11 3 BGEQ 46s 3$ : Check for high bit of C(n) set | 
08 A? D6 08 11 INCL 8(R7) ; Bit set. Adjust t C(n-1) 
FFEA 57 FFFFFFFC 8F 54 F1 08 1124 3$ ACBL R4, #-4, R7, 2$ ; Loop until C¢ 5) though C<(3) are in 
tt 1125 ; 5 though T8 
O83F 1126 
O83F 1127; 
O83F 1128 ; Generate the low 128 he - he product K*C'' = L. This product is 
O83F 1129 ; equivalent to multipl ying imes C'* modulo 8. The result of the 
O83F 1130 ; 135 it is in 14/17 with bits 31:29 of T4 the octant bits, and the ressining 
O83F 1131; its the faction bits of the product. The last 53 fraction bits (bits 
083F 1138 : 0 of T6 and 31:0 of T5) are non-valid tent ten bits that will be used 
O83F 1135 ; later if more fraction bits need to be generated. 
O83F 1134; 
O83F 1135; : : 
O83F 1136 ; Multiply the high order bits of K (RO) times C'' and store the result in 
O83F 1137 ; TO/T2. 
O83F 1138 
6E 00 04 A4 50 7A O83F 1139 EMUL RO, 4(R4), #0, (SP) TO/T1 = KHI*C(3) 
04 AE 04 AE 08 A4 50 7A 0845 1140 MUL RO, 8(R4), 4(SP), 4(SP) ; TO/T2 = KH EC C2) eC (3)] 
64 OC AS 50 C5 084D 1141 MULL3 RO. 12(R4S, (R4) ; 14 = Low 32 bits a KHI*C(1) 
08 AE 64 C0 0852 1148 ADDL2 (R4), 8(SP5 ; TO/T2 = KHI*C'* modulo 8 
0856 1144 ; Multiply the low order bits of K (R1) times C** and store the result in 
0856 1145 ; 1T4/T8. 
0856 1146 
00 04 A4 51 7A 0856 1147 EMUL R1, 4¢R4), #0, (R4) ; T4/TS = KLO*C(3) 
04 A4 04 A4 08 A4 51 7A Baas 1148 EMUL R1, 8(R4), 4(R4), 4(R4) : 14/76 = RLOSCECTD :€(2) 3 
08 AG «08 AS «(OC AG CSSA 0 64 1149 EMUL R1. 12(R45, 8(R4S, B(R4); T4/T7 = KLO®EC(1):C(2):C(3)) 
65 51 C4& O086C 1150 MULL R1, (RS) ; T8 = KLO*C(0) 
Oc Ad «6065S O86F 1131 ADDL2 (R5), 12(R4) : 14/17 = KLO®C'' modulo 8 
837 1138 3; Add KHI*C'* to KLO#C"' to get K*C'' Store the result in 14/T7 
04 A4 6E CO 087 1155 ADDL2 (SP), 4(R4) : 
08 AS 04 AE DB 087 1136 ADWC 4(SP$, 8(R4) : 
OC AS «08 AE d8 Baas 1188 ADWC B(SP), 12(R4) 3 T4/T7 = K*C'* modulo 8 
est 1159 ; 
a8 1 1160 ; At this point there may or may not be enough | valid bits in R3/R4 to wrryrete 
0881 1161; Y. If the first 12 fratction bits are all 1's or O's, there a possibility of 
1 1196 ; loss of slenificance when computing ¥ Consequently. we must check for loss 
! 1197 : of significance before converting 14/T7 to Y and I. | 
1 1165 ° 
65 FC AS 00008000 8F (C1 1 1166 ADDL3 #*X8000, -4(R5), (RS) : If the first 14 fraction bits are 1's | 


—- 
oo 
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eons ke : 


65  3FFFO000 ef 


vo 


-o° 
row 


BITL #*X3FFFOOOO, (R5) and the requced arg = 1-f or the 
BNEQ CONVERT first 14 bit are 0 and the reduced 


:16 are 0 and significance will 
be lost. 


; More bits need to be generated to cover the loss of gs Hee - + There are 
not enough registers to hold all the potential entra bits, so that the bits 
already generated must be put on the stack. 


QOOOO9CS*EF 16 JSB GEN_MORE BITS 


7 

| 

359 = f, then (and only then) bits 
9 fic | 
: Generate 85 additional bits and add | 
| 

| 


a P 
E E 
A 1167 
1198 
16 
1170 
1171 
1126 
117 
1174 
1175 
1126 
117 
gt 8 
0899 1180 ; them to existing bits. Results are 
0899 1181 ; . stored in 13/17 
54 046 C2 0899 11 ¢ SUBL2 #4, R4 ; Adjust R4 to reflect the addition of 
089C 113 : _ another longword of K*C'' 
8 #29, -4(R5), ; Check oss of significance is due 
15 FC AS 1D €O O89C 1184 BS 29, -4(R5), 4$ if | f signifi i 
: Hy 3; to leading ones or zeros | 
8 } 1185 ; Lost significance due to leading zeros | 
65 10 A4 OF 00 EA QBA1 1189 FFS #0, #15, 16(R4), (R5) |; If at least one bit is set. This 
1 12 QBA7 1190 BNEQ co 3 _ means lost significance was minor. 
OC AG OEFFFFFF 8F D1 OBAD 1191 CMPL M*XOEFFFFFF, 12(R4) ; If one of the three high bits is set, | 
BLEQ CONVERT $ ost significance was minor. 
0086 3 6 : 119 BRW LEADING_ ZEROS spe ; 
| 
6 1194 = 
; Lost significance due to leading ones 
o ? 1132 L ignifi d leadi 
65 10A4 OF 00 0886 1197 4$ FFC #0, #15, 16(R4), (R5) ; If at least one bit is clear. This | 
Oc 12 O8BC 1198 BNEQ CONV ;  meanslost significance was minor. 

OC AS F8000000 BF D1 O8BE 1199 CMPL #*xF8000000, 12(R4) : If one of the three high bits is | 
02 1€ O8C6 1500 BGEQU CONVERT ; clear, lost significance was minor. 
30.—=Cié«é#?M OBR \$8) BRB LEADING_ONES | 

OBEA 1303 
St ! Be CONVERT: 
Q8CA 1206 ; Isolate octant bits and convert fraction bits to a pair of D-format | 
CA quantities YHI and YLO 
B13 =n | 
65 FC AS QO3 1D €EF 6 CA 1209 EXTZ2V #29, #3, -4(R5), (R5) ; T8 = octant bits 
FC AS £0000000 ar CA 0 00 1210 BICL #*x£0006000, =4(R5) : Clear octant bits | 
54 cs ae D8 121 SUBL #12, R5, R4 7 R4 points to low order bits of h 
OOOO0A0B'EF 16 8D 1 Ie JSB CVT_T0_O0UBLE : RO/R1 = 2°29*h_Lo | 
E2 121 ; R6/R7 = 2°29eh_hi 
56 0100 8F A2 3 1214 SUBW2 #*X01D0, R6 3 R6/R7 = h_hi 
2 14 E7 1215 BGTR 3$ ; Check for h_hi =_0 
6 04 O8E9 1 18 CLRL 6 ; Restore h_hi to 0 
4 f bees : ' 3$ dh * ; Check for h_lo = 0 
50 0100 BF Ag a 15 SUBW2 PO a RO ; RO/R1 = h_lo . 
0965 =«€ FS : 1$: BLBC —_- (R5), 2 : Check forodd or even octant bits | 
; ctan s are odd. Reduced argument equals 1 - h. 
rf i Octant bit dd. Reduced t ls 1 h 


- 


XSGSINCOS ; Floatin 
eu b0e REDUCE_LA 
56 43FD OBF7 
8S S8 Oe 
9 
9 
52 20 AE 00 09 
OOAA 31 0904 
907 
907 
907 
907 
907 
907 
907 
0907 
0907 
907 
0907 
0907 
0907 
0907 
0907 
0907 
0907 
0907 
0907 
0907 
0907 
0907 
65  000003D0 8F DO 0907 
090€ 
090€ 
OC AG FFFFOOOO BF «D1 s«O90E 
1A 0916 
QOO009CS'EF 16 0918 
OC AG FFFFFFFF 8F D1 O9IE 
1E 1A 0926 
09 8 
08 AG 0406 7 09 
64 FC A4 7D 92D 
FFD3 65 0200 BF 0400 8F 3D 0931 
0938 
0938 
09 B 
9 
9 
938 
938 
338 
338 
938 
9 
938 
938 
53 04 0938 
52 OB AE O03 1D EF 0930 


) 
E 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
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80 


6-SEP-1984 


SUBG3 Re. #1, R6 
MNEGG RO, RO 

H Get octant bits 

2$: MOVL 32(SP), R2 
BRW GET_YH{_YLO 


t 
sign 
the 
bits to compute Y are available. 
allo 


to Y. 
LEADING_ONES: 


; If processing continues here it is known that the loss of significance is due 


; to a string of leading ones. 


MOVL #L_INT_WEIGHT, (RS) ; 
LOOP_1: CMPL #*XFFFFOOO0, 12(R4) ; 
BGTRU CONVERT_1 : 

SB N_MORE_BITS : 

CMPL #-17 12(R4) : 

BGTRU CONVERT_1 : 

MOVQ 4(R4), 8(R4) : 

MOVQ  =4(R4S, (R4) : 
ACBW #W_MAX_WEIGHT, #W_TERM_WEI 


; argument is either not 
n F-format (i.e. underf 


; argument to signal this occurance. 


; returned in Ri. 


CLRL B 
EXT;V #29, #3, B(SP), R2 


his point it has been determined that there is a major loss of 
ificance and the processing begins a looping phase. 
Loop will generate additional extra 


The wots one factor is groqter then 1024. 
distinguishable from 1 or too small 
low.) Zero is returned in R4 for the reduced 
Note that under these conditions the 

correct function value is one of the values 0, +/-1. The : 
; correct choice is determined by the calling program based on the octant bits 


AX/VMS Macro Vv04-00 
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nent =Y=1 = h hi 
RO/R1 = -h_lo 


Page 28 | 
. 48) | 

| 

| 

| 


R2 = octant bits 


Each iteration of 
bits of K*C’ until enough significant 


RB gy AR Ane go ore this time the nine longwords 
used as follows: 
TO/T2 Temporary storage used when generating extra bits. | 
13/17 Contains all significant bits generated so far. 
T8 Contains a counter, W, indicating the appropriate exponent 


of the last longword of fraction bits used in converting 


T8 = -~ bias for last longword 
of the product K*C' 


Check for enough significant bits 

Enough bits. Convert to floating. 

T2/T/ contains K*C*' 

Check for all 1's 

Not all 1's. Enough precision bits 

to compute Y 
| 
| 


Compress representation 
Ke ee 


) 

GHT, (R5), LOOP_1 

Increment weighting factor. If 
weighting factor is greater than 
1024 then no more bits need to be 
generated. 


This means that the reduced 
to be represented 


Reduced argument is zero 
Re = octant bits 


D.13 
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G 8: MTHRTL.S RCIU V¥SSINCO.MAR; 1 
BRw RESTORE 


CONVERT_1: 
ADD 


atin 
E_LA 
0078 31 43 
ie 
54 06 C0 ace L2 #4, RG ; oe a fe*ats.te ce ar bits of h 
QOOOOADB'EF 16 094 Cvf_TO_pOUBLE : 
94F : RORY = “uehch 
56 FEOS CF 46056 43FD «(094F SUBG3 6, 6 1° T0_32, R6 : RO/R7 = 24We(T = h_hid 
Q 50 52FD 0956 MNEGG RO, R ; RO/RI = =2*We 
52. FC AS 0 1D «EF «6OO95A EXTZV § o. =4(R5), R2 : R2 = octant bits 
6 65 C2 096 SUBL (R55, R6 + R6/R7 = 1 = 
50 ° 96 TSTwW ; Check for h_l6 = 0 
4A 96 BEQL GET_YHI_YLO : 
50 65 C2 096 SUBL fs : RO/R1 = = h_lo 
45 11 096A BRB GET_YHI_YLO 
096¢ 
96C 
096C LEADING_ZEROS: 


) 
E 
} 1 
1588 
1284 
} 5 
1289 
1 50 
1291 
1 5 
129 
1294 
13 | 
1599 
ge 1a 
096C 1300 : If processing continues here it is known that the loss of preryr cance is due 
096C 1301 ; to a string of leading zeros. Note that it is known that the loop for 
096c 1 oe : leading zeros will terminate before an underflow condition occurs so that the 
Boek 1307 ; loop does not include a test for underflow. 
65 00000300 8F DO 096C 1305 MO'YL #L_INT_WEIGHT, (R5) ; 18 = bias for last longword 
+44 1307 ; of t e product K*C' ¢ 
OC AS «QOOOFFFF 8F D1 0973 1308 LOOP_O: CMPL #*XOOOOFFFF, 12(R4) ; Check enough fraction bits G 
B 19 097B 1309 BLSS CONVERT_0O ; Enough bits. Convert to floating G 
000009C5 "EF 16 097D 1319 ace GEN MORE_BITS 3 T2/R contain K*C* G 
OC AS) =6—DS5 «60983 «1311 18: TSTL 12(R4) ; Check for all O's G 
10 12 0986 1 1¢ BNEQ CONVERT_O ; Not all 0's Enough precision bits. G_ 
08 AS «20046 AG) «67D SCO0988 S131 mMOva 4(R4), B(R4) ; Compress representat: on G_ 
64 FC AS 7D 098D 1314 MOVa  =4(R4S, (RG) : of Kec! G_ 
65 0200 8F AO 0991 1315 ADDW #W_TERM_WEIGHT, (R5) : Increnent weighting factor. G_ 
Te Ee hl Son § | is 
0998 1318 CONVERT -05 G 
54 046 CO 0998 1319 #4, R4 7 R4 gets s to low bits of h G 
OOOOOAOB'EF 16 44 : $9 188 Cvf_TO_DOUBLE : ares . sa ue hang | : 
52 FC AS O3 10 EF O9Al 1 é EXTZV #29, #3, -4(R5), R2 ; R2 = octant bits 
56 65) «6C2 «(O9A7 «1 SUBL2 (R55, R6 : R6/R7 = h_hi 
50 8B 9AA 1324 TSTW RO : Check for h_lo = 0 
0321 9AC 1325 BEQL GET YHI_YLO 3 
50 65 (C2 Q9AE : § SUBL2 (R5J, RO ; RO/R1 = h_lo | 
1308 GET_YHI_YLO: Pe = | 
1 mova R6, R4 ; R4/RS = high bits of 
1330 CLRL R?7 ; R6/R7 = high 24 bits Y = YHI 
1331 SUBG RO, R4 5 
1 ¢ ADDG RO, R4 ; R4/R5 = YLO 
1333 RESTORE: 
1334 Bisse 36(SP) ; Restore IV bit and exit 
: 5 oeee #40, SP ; Remove mask and temporary storage 
1339 | 
| 


E 13 | 
XSGSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 :05:25 VAX/VMS Macro Vv04-00 Page 
petit REDUCE. LARGE ae 94:93 6 [MTHRTL.SRCJUVXGSINCO.MAR:1 9 a 
9c5 (1 
9c5 1 § GEN_MORE _BITS: 
gH | 
9¢c5 «(1 tg 3 : This gubroys ine generates 85 extra fraction bits and puts them to the 
905 134 ; existing bits. This routine is always entered hg a JSB instruction. 
ae | te 3 3 sheet pmainade SP + LS to the first longword BEFORE TO, rather than T 
ge ie 
3c @ C2 43 ! rf SUBL2 #4, R2 ; Ad ron Fipubat a Ov581 4 next quadword 
56 620«COS3s79 3c 1350 ASHQ = RB, (RZ), Rb : R7 = 
17 =—*18 Bare ! 23 BGEQ 1$ ; eneh " high bit is clear 
O3te ! 35 ; Logic to process unsigned values greater than 2°31 - 1 
04 AE 00 51 7A 9CE 1355 EMUL R1, 87, #0, 4(SP) 3 
ag 51 co Bap8 1 2$ ADDL2 Ri, : TO/T1 = KLO*C(n) 
08 AE O8 AE 5 50 7A 0908 135 EMUL RO, R7, 8(SP), 8(SP) 3 
OC AE 50 COQ O9DF 1358 ADDL2 RO, 12(SP) 3; TO/T2 = K*C(n) 
OD 11 ead 1323 BRB 2$ 
945 1361 3; Logic to process unsigned values less than 2°31 
04 AE 00 57 51 7A Q9E5 1368 1$: EMUL R7, #0 (SP) 3 TO/T1 = KLO*C(n) 
08 AE OB AE 57 50 7A O9E 1364 EMUL RO: R7., BCEP). 8(SP) > T0/T2 = K*C(n) | 
aoe : 96 3; Add new bits to old 
04 Re OC AE OS O96 1360 AbwE’ ToRSPS, “ueRaD 
08 1E€ O9FB 1359 BCC i$ : Check for carry from previous add 
08 AS D6 O9FD 1371 INCL 8(R4) ; Propagate carry 
03 1€ OA00 1376 BCC 3$ 3; Check for carry from previous add 
OC AS D6 OAD2 137 INCL 12(R4) ; Propagate carry 
FC AS «604 AE DO OQADS 1374 3$ MOVL 4(SP), -4(R4) 3; Move new low order bits to end of 
OAOA 137? : of old low order bits 
05 OAOA 1376 RSB 3 
QAOB 1377 
QOAOB 1 8 
OAOB 137 | 
AS ia 
AOB 1 ¢ CVT_TO_DOUBLE: 
AE | 
QA0OB 1385 : This routine converts an array of three longword pointed to by R4 to a pair | 
AOB 1 g: 3; of G-format values. The results are returned in RO/R1 (low 48 bits) and 
hoe : 3 gE RT (high 48 bits). ; 
50 84 sefo A0B 1 89 CVTLG (Rd), RO ; RO/R1 = Low 32 bits of h 
07 1a OAtl 1391 Berk : Adjust for signed | 
4 D6 OA! 1 3 INCL (R4) 3 conversion error 
63 1E QAIS 139 BCC 1$ 3; If necessary, 
04 AS D6 OAI7? 1394 INCL 4(R4) : «propagate carry 


F 13 
XSGSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-19 AX/VMS Macro v04-00 Page 31. 
oboe REDUCE. LARGE gre ee 1 9be 9Fi8a:8 PATANte. caedovNcsinco.man:1 2% (ip. 
A2 OAIA 1395 1$: SUBW #W_TERM_WEIGHT, RO RO/R1 = (low 32 bits of h)/2*32 | 
FC AG 64 +P FOSS ef £8 AIF 1 38 2$: sict§ PRET POODD. (R4), <4 (R45 
$4 FC AG C2 «(0A28 «139 SUBL2 =4(R4), (R4S : 
e. oS AG ri 4 Aec 1 38 CVTL -4(R4), R2 3 
g 4OF A31 139 ADOG Re RO 3; RO/R1 = (low 48 bits of h)/2*32 
1 A ; 1400 BEQL 3 3 
50 0 F A A 1401 SUBW2 #W_TERM_WEIGHT, RO ; RO/R1 = (low 48 bits of h)/2%64 
4 “ere A 14 4 3$ CVTLG (RZ)+, R2 ; R2/R3 = next 16 bits of h 
25 AG 14 BEQL 5$ 3 
02 16 QOAS2 1404 BGTR ; Ad ues for signed conversion error. 
64 06 AS& 1405 INCL (R4) 3 e that no carry is Rasy, sible 
52 200 8F A nate 1406 4$ SUBW2 #W_TERM_ WEIGHT, R2 3 R2/R re (next 16 ots of h)/2*32 
6 646 GEFD AGB 14 5$ CVTLG (RZ), RB ; R6/R7 = high 32 bits of h 
06 18 OASF 1408 BGEQ 3 wae for signed conversion 
56 FD64 CF 4OFD OAS1 1409 ADDG G_2_10_32, R6 3 
56 52 40FD OAS7 1410 6$ ADDG R2, R6 3 ROIR? = = (high 48 bits of h)/2*32 
OF OAse late oie 


> + aaaa 
-004 


es Ww wn" 
o 
ouMl MOOW® woOwdw wWwodwW 


“ 
— 
Sa] 
m” 
aa) 
~ 
MOVWM WEN 
OW WOOF WwWTO"SO ONTOS 


wu 
o 


“ 

™~ 
Ww 

w 

~™ 
Ww 
mur 

~ 

Oo 


Cad 
mo 
mo 


eee 
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POWER_MOD_360_0: 
.QOR 


ee ee ee el el ee cel ee el el eel eel el eel eel el eel el cece eel el eel eel el cel el cel el el el el el el cel cel eel eel cel eel eel eel el eel el eel cel el el eel eel el el eel 
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G 13 
nt Sine, Cosine and Sincos at det 


“SEP=19 
-SBTTL REDUCE _DEGREES 


This routine assumes that the absolute value of the argument is in RO/R1. 
The reduction process is performed in two 
the reduction reduces the argument modulo 
and the second stage reduces the argument modulo 45 to a value less chen 45. 


; Constants used in this reduction: 


ORD - * 

» WORD 16, ; ‘ 64; 

«WORD 256, 152. 304, 
POWER_MOD_360_1: 

.GORD” 136, 573° 

~ WORD 16, 35° 64, 

WORD 256, 152, 304; 
REDUCE_DEGREES: 

CMPW #*x4360, RO 


BGTR —_LAST_STEP 


AX/VMS Macro V04-00 


Pp 
MTHRTL.S REI UVXGSINCO.MAR; 1 


tages. The first stage 
60 to a value less that 


; Powers of 2 modulo 360 for t1 = 0 


Powers of 2 modulo 360 for t1 <> 0 


Compare ix! with 2°53 
Branch to special logic for med arg 


It is assumed here that the argument is greater than 2°53. 


; The ar maaan is reduced as follows: 


J= 64 - Now 
(2*12)*2 = (2*9)#(2*15) = (2 “9) « ( 
= 2" aseist2)- 


not zero, K = 2* 
For tl =0K =2 
and congruent to 


K'#(64*J1 + J2) is congruent to s 


OVL RO, R2 
BICL2 #*X7FFO, RO 
1SW #*X4350, RO 
SUBL RO, 
MOVQ 


RO, R3 
BICL @*XFFFFSFFF, R1 


SUBG RO, R 

suBu3 dt of RO 
ADDG R3, RO 

ROTL #-4, ; a 
DIVW3 He: 

MULW 

SUBW 

TSTW R3 

BNEQ 1$ 


; Check for t1 
accordingly 


< f.< 1. And let J = 2*53*f = 


/ 

pince, 2°30 = 64 nodyle 360, we have that 
= i, - Po aesee + t2. Note that 
$5 12 mod 360. Hence, if t1 i 
. a ere )«( 2negy! : iSos2*¢5 agduto zoo. 
TF Consequently, define Kk’ congruent to 2*t2 if tl 
156*2*t2 otherwise, where 0 =< K ‘ 

modulo 360 and x’ me 2°53. 


Then x’ 


high longword of X 
. exp bits of X 


nunnnw 
= 
™m 


J1 
gar) + J2 = J modulo 45 


ce 
—_ «= 


ro— 
nor 
a 
oo 
_ 


= 0 and choose K' 


] 
H 13 
XSGSINCOS ; Floatin nt Sine, Cosine and Sincos 16-SEP-1984 :05:25 VAX/VMS Macro V04-00 Pa 3 | 
petty kepuce DEGR aes 95:93:82 MTHRTL. SRCJUVXGSINCO.MAR;1 9” 30) | 
52 88 AF42 4DFD OACF CVTWG POWER_MOD_360_0CR2], R2 : R2/R3 = K' 
0 11 M03 BRB $ ; 
52 28 AF42 4DFD OAD 1$: CvTW Spee MOD_360_1CR2], R2 : R2/R3 = 
5 52 44FD nee 2$: MULG R2, RO 3; RO/R1 = she (mod 45) 0 =< RO < 2°55 
AE1 
nt LAST _STEP: 
wt Argument reduction scheme for arguments with absolute value less than 2°55 
QAE1 : The reduced ,orqueent Y is computed as follows: 
AE1 3 Let 1 = 
AE1 $ if I 2 pe. 
AE1 3 then Y = X = 45e] 
OAE1 3 else Y = tis1}245 - x 
Oar] 
50 42460 8F B81 OAE1 CMPW #*X4240, RO ; Compare 7°36 with iX! 
8 18 OAE BGEQ NO 1CVEREL OW : 
56 50 F563 CF 45FD OAE MULG3 G_T_OV_45, RO, R6 3 R6/R7 = {X1/45 


: Turn off IV to avoid an exception in EMODD 


RIPINOPONININS 9 9 SS “DOO OOOOOCOO COO OOO OOO OWWWWWOWOOWOWO NNN We: 


AFAR “OOODNOUSWN “OOD NAUES WN $9 ODONAU EWN 0 OD NAME WN O OONOUS wn 


52 DC OAEF MOVPSL R2 ; Move PSL to R2 
52 FFFFFFDF 8F CA OAF1 BICL #*C<PSLSM_IV>, R2 ; Save current IV bit 
20 are BICPSW #PSLSM_IV ; Turn off integer overflow trap 
54 53 56 00 08 54FD OAFA EMODG #1, #0, R6, R3, R4 ; R3 = low 32 integer bits of ixXi/45 
49 3; R4/R5 = fractional part of 'X:/45 
52 B68 pest BISPSW R2 ; Restore IV bit 
56 54 42FD 0803 SUBG2 R4, R6 ; R6/R7 = Integer part of 'X:/45 = 
2F 53—OED 07 BLBC RS, EVEN ; 
56 08 40FD OBOA ADDG2 #1, R6 ; R6/R7 = 1 + 1 
16 11 443 BRB oD 
081 
0810 
0810 NO_OVERFLOW: 
53 50 Cl6C BF F538 cr 54FD 0810 EMODG G_1_0V_45, #X_1_0V_45, RO, R3, R4 
0818 ; R3 = I = integer part of !X!/45 
1753 =—«E9 dt LBC R3, CVT : Branch it octant bits are even 
56 53 9 C1 ig ADDL3 #1, R3, R6 ; R6 = 1 + 1 
56 «456 GEFD 08 CVTLG 6, RO > RO/R7 = he | 
56 F50D fF 44FD 6 ODD: MULG G 45 R6 3: R6/R7 = 45*(141) 
56, , 30 42kD 089 SUBG2 RO, : RO/R | 
53 8F aA 1CB #°XFB, RB ; Save onke last three octant bits 
5 4 RSB | 
53 4EFD : CVT: CVTLG R3, R6 : R6/R7 = 1 
gf 44FD 9 EVEN: MULG G n45 R6 > R6/R7 = -45*] 
0 40FD F ADDG2 RO, Rb > RO/R7 = 
F8 8F BA 43 BICB #°XFB, R3 : Save only last three octant bits 
05 re RSB 


yegesincos ielMERSSENEDITE Stee Conte and'sicon e-s6e-1Hs F7G5:35 ARAYA MBER anes 29" chy 


0848 1527 


J.13 
X$GSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 02:05:25 VAX/VMS Macro v04-00 Page 35 | 
S006 REDUCE. DEGREES SrSep=1986 11:98:86 MTHRTL.SRCJUVXGSINCO.MAR;1 (By 
48 1529 | 
6 : 9 -SBTTL RADIAN_POLYNOMIALS : Polynomials for arguments in radians | 
48 1 
ais 1538 | 
48 1535 ; Polynomial evaluation for DCOS(Y) for Y in radians 
te 1889 ° | 
48 1338 P_cos_R: 
53 36 8000 8F AB 48 1 BICWw3 #*x8000, R6, R3 3 
5 4000 ef } pB% 1309 CMPW #*x4000, R ; Compare 1/2 with {YHI! 
6 4 53 154 BGTR LEQL_HALF ; Sufficent overhang is available 
085 1308 NEEDS_DOUBLE: 
7E 2 7D 085 154 MOVQ R4, -(SP) 3; Save YLO 
7E 54 6 41FD OBSB 1544 ADDG R6, R4, -(SP) 3; Save Y 
50 6e gE 45FD O0B5D 1545 MULG3 (SP), (SP), RO : RO/R1 = Y 
F607 CF 0 56 55FD pBbe 1368 POLY RO, #COSLENR2-1, COSTBR2; RO/R1 = Q(Y*2) 
54 56 BE 41FD 0869 154 ADDG3 (SP)+, R6, R4& : R4/RS5 = Y + YHI 
54 of bas | pass 1548 MULG2 (SP)+, R4 3; R4/R5 = YLO*(Y + YHI) = A2 
0 1 OB72 1549 BEQL 3; Check for A2 = 
54 10 A2 0B74 1550 SUBW #*x0010, R4 : R4/R5 = A2/ 
a8 ae bara Opre ieee ure RE. Re ape mies fee 
56 10 = A2s(OBF 1388 SUBW2  #*X0010, RO + R6/R7 = YHI*2/2 
56 ©6008 ry 4 OB82 1554 SUBG #1, R6 3; R6/R7 = -(1 = YHI*2/2) 
50 56 42FD 0B86 1555 SUBG R6, RO ; RO/R1 = GCOS(Y) 
0 Opee 1889 ~ 
0B8B 1558 LEQL_HALF: 
se secu ABP tp aes ae a MOA Ya 
F596 CF 07 56 S5FD ti 1561 POLYG R6, #COSLENR1-1, COSTBR1; RO/R1 = GCOS(Y) 
05 OB9A 1366 RSB 
0898 156 
Be ie 
0898 1298 3; Polynomial evaluation for -GCOS(Y) 
98 1567; 
38 1369 N_COS_R 
53 56 8000 8F AB 98 1390 “~~ Bicw3 #*x8000, R6, R3 ; ; 
53 4000 8F 61 Al 1571 CMPW #*x4000, R3 ; Compare 1/2 with ‘YHIi 
7 14 : § 1328 BGTR 2s 3; Sufficent overhang is available 
7E 4 7D ; 157 MOVQ R4, -(SP) ; Save YLO 
43 54 56 41FD BAB «(1574 ADDG R6, R4, -(SP) ; Save Y 
50 6 : 45FD OBBO 1575 MULG3 (SP), (SP), RO + RO/R1 = 
F5B4 CF 2 55FD 1376 POLY RO ACOSLENR2-1, COSTBR2; RO/R1 = Q(Y*2) 
54 56 BE 41FD 157 ADDG3 (SP)+, R6, R4 > R4/RS = Y + YHI 
54 “sre 1 4 MULG (SP)+, R4 : R4/RS = YLO®(Y + YHI) = A2 
0 1 157 BEQL 1 3; Check for A2 = 0 
4 1 A2 1580 SUBW #*x0010, R4 : R4/R5 = A2/ 
0 3 42FD 1 1 18: SUBG R4, RO ; RO/R1 = Q(Y¥*2) = A2/2 | 
6 _ i § MULG R6, R6 : R6/R7 = - 
6 10 A 1 SUBWo © #*K0010, R6 ; RO/RZ = YHI*2/2 | 
28 9 t¢fe 1584 SUBG #1, R : R6/R7 = =-(1 = YHI*2/2) 
50 6 45FD 1585 SUBG RO, R6, RO ; RO/R1 = -GCOS(Y) 
| 


“ole v04-00 Page 
i =SEP=1984 02:05:25 VAX/VMS Macro 8:1 (38) 
: Point Sine, Cosine and Sincos 16-SEP- . 17) RTT SRE SOvacsIocO.MAR: 
So pene habLan POLYNOMIALS 3 bolynomials for arg 6=SEP 1382 %¢ F 
; oF s : R6/R7 = Y 
DF 1 2$ ADDG R4, R6 : AOUR? = Y*2 
se é garb ES ; POL YG Re. ACOSLENRT~1, COSTBRI; nora : GCOS(Y) 
ino gf g ; ne E 1591 XORW #*x8000, RO : R 
. repeal Os 4 9 RSB 
8 Bs 
4 He Polynomial evaluation for -GSIN(Y) 
F4 138 . 
fo] 35 oe #*xX8000, R4 : 
000 8F AC F4 159 XORW * - ; aise? w <4 
4 000 8F AC bare 1601 XORW #*x8000 
Oar 1808 Polynomial evaluation for GSIN(Y) 
OSE 1805 ieee 
i, -(SP) 3; Save 
re Ht yf SE Tet nova, ae, 2450), : fave 1 
7E 56 54 41FD OCO 1608 eS) tsp) Ue aes f 
ef 65 gf eefe tats 1303 poty oes #SINLENR=1, SINTBR ; pd . ht 
FS9E CF 06 Re, ASIN ; ne 
PREP EL | Ba Ghee | RSLS 3 Bie carers 
35 SE tOED be : ADDG2 RO, RO : 
50 56 40FD OCIA 161 a 
BE ek 
1F 1616 
Sit 1617 
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Sine, Cosine and Sincos 16-SEP-1984 
S ; Polynomials for argu 6-SEP-1984 
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+ Polynomial evaluation for GCOS(Y) for Y 


G_°-0V_PI. R6 


R4, =(SP) 

R6, R4, =(SP) 

(SP), (SP), RO 

RO, MCOSLENC2-1, COSTBC2 
(SP)+ RG 


’ ° 


i R4 
#*X0020, R4 
R4, RO 


R6, R6 
#*x0020, R6 
#1, R6 
R6, RO 


R4, se ; 
R6, WCOSLENCI~1, COSTBC1: 
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», (SP), RO : 
RO, #COSLENC2-1, COSTBC2: 
(SP)+, R6, R4 : 
(SP)+, R4 

#7 x0020, R4 

R4, R 

R6, RS 

#°X0020, R6 

RO, R6, RO 


R4, R6 
R6, R6 
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in cycles 


Compare 2/pi with !YHI! 
Sufficent overhang is available 


Save YLO 
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; Compare 2/pi with ‘YHI: | 
Sufficent overhang is available 


Save Y 
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AX/VMS Macro V04-00 
MTHRTL.SRCJUVXGSINCO.MAR; 1 


Polynomials for arguments in cycles 


Y + YHI 
YLO*(Y + YHI) 
2=0 


A2/4 
Q(Y¥*2) = A2/4 
YHI“ 


Y + YHI 
YLO#(Y + YHI) 
or A2 = 0 


A2/4 
Q(v*2) = A2/4 
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VXSGSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 :05:25 VAX/VMS Macro v04-00 Page 
BuO0e EYCLE POLYNOMIALS ; Polynomials for argu ~$Fb= 1382 99:93:88 MTMRTL.SRCJUVXGSINCO.MAR;1 © 38) 
F529 CF 7 SSFD CB8 «167 POLY R6, MCOSLENC1-1, COSTBC1; RO/R1 = GCOS(Y) = 1 
50 00000000 86006010 +t 28 a3r cet 187? Cun RB” #-1, RO 5 AA = -GCOS(Y) 
ccD 1878 
CCD 1680 ; 
CCD 1681 ; Polynomial evaluation for -SIN(Y) 
CCD 16 § 3 
CCD 16 
CCD 1684 N_SIN_C: 
54 8000 8F AC OCCD 1685 XORW #*x8000, 3 
56 8000 8F AC De 16 § XORW #°x8000, : RG/R7 = - Y 
CD7 i 3 : 
CD? 1689 ; Polynomial evaluation for GSIN(Y) 
OCD? 1690 ; 
OCD7 1691 
OCD7 1698 P_SIN_ 
7E 54 7D OCD7 169 R4, - 3; Save YLO 
7E 56 54 41FD OCDA 1694 ADDG R4, R6, al ; Save Y 
54 6€ 6€ 45FD OCDF 1695 MULG3 (SP), (SP), R4 : R4 = Y* 
FS7D CF 06 54 5SFD OCE4 1696 POLY R4 #SINLENC™ -1, SINTBC ; RO/R1 = P(Y*2) 
50 BE 44FD OCEB 1697 MULG2 (SP)+, : RO/R1 = Y*P(Y*2) 
50 6E 40FD OCEF 1698 ADDG (353 an ; RO/R1 = YeP(Y*2) + YLO 
6E 20 A2 OCF3 1699 SUBW #10620, (SP) : (SP) = YLO/4 
50 E 42FD OCF6 1700 SUBG (SP) RO ; RO/R1 = YeP(Y*2) + 3/4*YLO 
54 56 7D OCFA 1701 MOVQ R6, : R4/RS = YHI 
54 20 A2 OCFD 1786 SUBW #40020, R4 3; R4/RS = YHI/4 
56 54 42FD O0D00 170 SUBG R4, R6 : R6/R7 = 3/4*YHI 
50 56 40FD OD04 1704 ADDG ‘ ; RO/R1 = GSINC(Y) 
05 0D08 1039 SB 


2 ad 
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N 13 
VX$GSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 02:05:25 VAX/VMS Macro V04-00 Page 39 — 
Petty bEGREE. POLYNOMIALS o- See 138 94:93:83 MTHRTL.S REI VXGSINCO.MAR; 1 . Bay | 
D09 1708 .SBTTL DEGREE_POLYNOMIALS | 
D09 «1709 | 
D09 1710 
DO9 1711 P_COS_D: | 
56 F352 CF S1FD 0D09 216 CMPG 90_OV_PI, R6 ; Compare 90/pi with Y | 
40 18 OQDOF 171 BGEQ g ; Double preg 'sion isn't needed 
54 56 56 45FD OD11 1714 MULG3 = R6, + RO/R1 = 
F5C3 CF 07) 54 S5FD «0D16 «1715 POL YG R4, ; CSstENT. COSDTB1 : RO/R1 = Q(Y*2) 
54 56 00070000 BF cB Dip 1216 BICL3 #*X70000, R6, R4 : 
55 04 0025 1718 CLRL_ RS : RG/RS = YHI 
52 36 34 43FD OD2? 171 SUBG3 4, R6, R2 : R2 = 
6 54 40FD ODec 1720 ADDG¢ R4, R6 + R6/R7 = Y + YHI 
56 52 44rd 0D30 1721 MULG2 R2, R6 : R6/R7 = YLO®(Y + YHI) = A2 
05 13 0034 17 : BEQL : Check for A2 = 
56 0000 BF A2 9036 17 SUBW2 #*Xx00D0, R6 + R6/R7 = A2/2°13 
50 56 42FD O0D3B 1724 1$ SUBG2 R6, RO + RO/R1 = Q(Y*2) = A2/2*13 
54 54 44FD «(OD3F «01725 MULG2 R4, RG + R4/RS = YHI* 
54 0000 BF A2 OD43 1726 SUBW #*X00D0, R4 3 R4/RS = YHI*2/2°13 
54 «(08 43FD 0048 1727 SUBG2 #1, R4 + RG/RS = -(1 = YHI*2/2°13) 
50 «54 42FD 4OD4C «(1728 SUBG R4, RO ; RO/R1 = GCOS(Y) | 
05 0D50 1729 RSB 
0051 1730 | 
56 56 44FD 0051 1731 2$ MULG2 6, R6 > R6/R7 = Y*2 
08 13 0D55 1732 BEQL : Check for Y = 0 
F542 CF 07 )«=656 SSFD «(0D57 ~«(17 POLYG R6, #COSDLN2, COSDTB2 ; RO/R1 = Q(Y*2) 
05 ODSE 1734 RSB 
ODSF 1735 
50 08 SOFD ODSF 1736 3$ MOVG #1, RO : RO/R1 = GCOS(Y) | 
05 0063 1737 RSB | 
0064 1738 
0064 1739 
0064 1740 N_COS_D: 
56 F2F7 CF SIFD OD64 1741 CMPG G_90_OV_PI, R6 ; Compare 90/pi with ¥ 
41 18 OQD6A 1708 BGEQ 23 ; Double precision isn't needed 
54 56 56 45FD OD6C 174 MULG3 6, R6, R4 + RO/R1 = Y*2 
F568 CF 07 54 55FD O0D71 1744 POLY Ra PACOSDLNI . COSDTB1 =; RO/R1 = Q(Y*2) 
54 56 00070000 8F CB OD78 1745 BICL x70000, R4 ; 
55 D4 OD80 1746 CLRL RS. > R4/RS = YHI 
52 56 54 43FD On8e 1747 SUBG3 R4, RO, R2 : R2 = YLO 
56 «54 4OFD «0D87 «1748 ADDG2 R4, RO + R6/R7 = Y + YHI 
5652 44rD 0D8B 1749 MULG2 R2, R6 : R6/R7 = YLO#(Y + YHI) = A2 
0S 13 OD8F 1750 BEQL : Check for A2 = 
56 0000 af A2 09) 1751 SUBW #x0000, R6 > R6/R7 = A2/2* 
0 56 42FD 0096 1738 1$ SUBG2 6, R : RO/R1 = Q(Y*2) = A2/2°13 
54 54 44FD ~=OD9A «175 MULG2 R4, R4 > R4/RS = YHI* 
54 0000 BF A2 oo9 1754 SUBWo #*KOODO, R4 > R4/R5 = YHI*2/2°13 
54 08 43rD DAS 1755 SUBG : : R4/RS = -(1 = YHI42/2°13) 
50 54 50 43FD ODA7 1796 $U853 RO, R4, RO + RO/R1 = -GCO 
aa a 
56 56 44ro DAD 1759 2$ MULG2 6, R6 : R6/R7 = Y*2 
Op 1 DB1 1760 BEQL  3$ : Check for Y = 0 
F4E6 cf 97 6 55FD ODB3 1761 POLYG RG #COSDLN2, COSDTB2 ; RO/R1 = GCOSD(Y) 
0 8000 8F aC DBA 126¢ XORW *x8000, RO + RO/R1 = -GCOSD(Y) 
ooh ie | 


B 14 
XSGSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 :05:25 YVAX/VMS Macro v04-00 Page 40 
yvase beGREE POLYNOMIALS ant et 9 94:93:83 ERTHRTL. SRE JUVXGSINCO.MAR: 1 . (24 
50 00000000 0000C010 8F 50fp pco 1763 3$: MOVG #-1, RO ; RO/R1 = GCOS(Y) 
0 pcc 1 § RSB 
pcdD 1 
ocd 17 $ N_SIN_D: 
56 56 S52FD 45 178 ‘1 MNEGG R6, R6 3; R6/R7 = -¥ 
50 56 56 4568 DDI 1559 a “ MULG3 R6, R6, RO ; RO/R1 = ¥*2 
4} dds 1228 BEQL RETURN 
F541 CF 6 0 SSFD DD 177 POLYG RO, #SINDLN, SINDTB ; RO/R1 = PCY*2) 
0 6 44FD ODDF 1774 MULG R6, ; RO/R1 = YeP(Y¥*2) 
56 060 BF Ae 434 1775 SUBW #*x0060, R6 3 R6/R7 = Y/2* 
0 6 4OF DE 1006 ADD R6, RO 3; RO/R1 = GSIN(Y) 
0 DEC 1777 RETURN: RSB 
ODED 1778 


C 14 
XSGSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 02:05:25 VAX/VMS Macro Vv04-00 Page 41 
Patty BEGREE. POLYNOMIALS g-8 P=1984 94:93:82 MTHRTL.SRCJUVXGSINCO.MAR; 1 ° (26) | 
DED 1780 | 
BED 178) 
444 1f 5 -SBTTL DEGENERATE_SOLUTIONS 
DED 1784 P_ONE: 
50 08 50fe DED 1785 MOVG #1, RO ; Answer is 1 
0 DF1 17 § RSB 
DF 17 
DF2 17 : 
DF2 1789 N_ONE: 
50 00000000 0000C010 8F sore DF2 1790 MOVG #-1, RO ; Answer is -1 
0 DFE 1791 RSB 
DFF 1036 
DFF 179 | 
DFF 1794 UNFL: 
re 1795 ; 
Bore 1736 ; Underflow; if user has FU set, signal error. Always return 0.0 
52 DC ODFF 1798 MOVPSL ne ; R2 = user's or jacket routine'’s PSL | 
00000000 ° GF 9° FB QE01 1799 CALLS #0, G*MTHSSJACKET_TST ; RO = TRUE if JSB from jacket routine 
04 50 £9 0£08 1800 BLBC ; branch if_user did JSB 
52 04 AD 3C OE0B 1801 MOVZWL SFSW_SAVE_PSW(FP), R2 ; get user PSL saved by CALL 
50 D4 OEOF 1008 10$: CLRL RO ; RO = result. LIBSSIGNAL will save in | 
0E11 180 ; CHFSL_MCH_RO/R1 so any handler can 
0E11 1308 : = fixup 
0D 52 06 €1 OE11 1805 BBC #6, R2, 208 ; has user enabled floating under f Low? 
DD OE15 1396 PUSHL P ; yes, return PC from special routine 
7E O0O°8F 9A OE17 180 MOVZBL #MTHSK_FLOUNDMAT, -(SP) ; trap code for hardware floatin | 
QE1B 1808 ; underflow convert to MTH$_FLOUNDMAT 
Q0E1B 1809 ; (32-bit VAX-11 engeptten code) 
00000000'GF 02 FB OE1 1319 CALLS #2, G*MTHSSSIGNAL ; signal (condition, PC 
05 OQE22 1811 208: RSB 3; return 
OE2 1316 
OE23 181 . END 
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= Wns 
oo t+#+epettt B6-O000 
22 aaa 


Dd 
Dd 
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Dd 
D 
1) 
$_ DOUBLE 


SIND 
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” +20 
So j= 
vaooow 
222222222 
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MOD 
MOD 
1V 
D 
L 
Ki 
$u_SAVE_PSW 
cos 


Ss 
S 
$ 
E 
N 
N 
N 
C 
C 
C 
STOR 
TURN 
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erereaeee 
gerernree 


; Floating Point Sine, 
gereraere 


COS 
“SIN 
HND 
TST 


were) 
nw 
<< 
I) 


ENERATE_ 
“ZEROS 
CF 
INT_WEIGHT 
HSSJACKET_ 
MTHSSJACKET_ 
MTHS$$SIGNAL 


L 
L 
" 


CVT_TO_ DOUBLE 


DEGENERATE 


LEADING 
LEQL_HA 
LONG 
LOOP 
LOOP 
L_co 
SIN 
T 


UVXSGSINCOS 
Symbol table 


DEG 


UVXSGSINCOS 
Symbol table 


SINTBC 
NTBR 


W_ADJUST 
WIMAX WEIGHT 
W"TERM_WEIGHT 
X_1_0V_45 


PSECT name 


. ABS . 
$ABSS 
_MTHSCODE 


Pass 


39338 bytes 


Initialization 
Command processing 


Pass 
ag table sort 


Symbol table output 
Psect synopsis output 
Cross-reference output 
Assembler run totals 


The working se 


(7 


Macro Library name 


-$255$DUA28: CSYSLIBISTARLET.MLB;2 
131 GETS were required to define 5 macros. 
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; Floating Point Sine, Cosine and Sincos i ak i dat 9 9¢:9 


GOOCCOCOOCOCOO 


toa mo san eases nases > 
H Psect synopsis ! 
Allocation PSECT No. Attributes 
00000000 ( 0.) 00 ¢ 0.) NOPIC  USR 
00000000 ( 0.) O1¢ +41.) NOPIC USR 
OOOO0E2S ( 3619.) 02 ¢ 2.) PIC USR 
or 
: Performance indicators ! 
ults CPU Time Elapsed Time 
33 00:00:00.08 00:00:01.14 
119 00:00:00.64 00:00:03.50 
208 00:00:05.88 OR: BB: 1B -83 
0 ba Bo 0:00:00.52 
318 00:00:03.91 00:00:10.83 
16 00:00:00.14 00:00:00.64 
3 00:00:00.02 00:00:00.02 
0 00:00:00.00 Be: 0a 8e-B2 
699 00:00:10.95 00:00:35.56 


Limit was 1650 pages. 
pages) of virtual memory were used to buffer the intermediate 
There were 20 pages of symbol table space allocated to hold 195 non-local and 
1873 source Lines were read in Pass 

10 pages of virtual memory were used to define 


+ 
' 
+ 


+ producin 


Macro Library statistics ; 


macros. 


Macros defined 


6-SEP=19 


CON 
CON 
CON 


6 object records in Pass 2. 
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ABS 
ABS 
REL 


8:86 


LCL 
LCL 
LCL 


AX/VMS Mac 
MTHRTL.SRC 


NOSHR NOEXE NORD NO 
R EXE RD 


NOSH 
SHR 


i symbols. 


EXE 


ro V04-00 
JUVXGSINCO.MAR; 1 
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WRT NOVEC BYTE 
WRT NOVEC BYTE 


RD NOWRT NOVEC LONG 
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BLSEYRSRES um seattstics + FUsetton Potet Sine Conion and Siicon SEES FREE HOTT MSSOVERG Mla muns 9%" cf 


There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESS1ON/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:UVXGSINCO/OBJ=OBJ$:UVXGSINCO MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) #MS 


i< 


2 


| 0265 AH-BT13A-SE ) NT CORPORATION 
gave VAX/VMS V4.0 D PROPRIETARY 


